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ABSTRACT 

A model is presented for the quasiperiodic component of magnetar emission during the tail phase 
of giant flares. The model invokes modulation of the particle number density in the magnetosphere. 
The magnetospheric currents are modulated by torsional motion of the surface and we calculate that 
the amplitude of neutron star surface oscillation should be ~ 1% of the NS radius in order to produce 
the observed features in the power spectrum. Using an axisymmetric analytical model for structure of 
the magnetosphere of an oscillating NS , we calculate the angular distribution of the optical depth to 
the resonant Compton scattering. The anisotropy of the optical depth may be why QPO are observed 
only at particular rotational phases. 

Subject headings: stars:magnetic field - starsmeutron - stars:oscillations 



1. INTRODUCTION 



iDuncanI (|1998f l proposed that a magne tar undergoing a 
iant fl are would b e set into oscillations. iTimo khin ct al.' 
200Cl| ): PTimokhinl ((2007) proposed that torsional oscilla- 
tions could modulate the local Goldreich- Julian charge 
density and noted that changes in the Goldreich-Julian 
char ge density are prop ortio nal to the harmonic numbers 
I, m. Ilsrael et"all (|2005D and lStrohmaver fc Watti (|2005l ) 
have reported detections of quasi-periodical features in 
the power spectra of SGR 1806-20 and 1900+40, which 
they interpret to be n=0 high-Z torsional oscillations of 
the crust. These features in the power spectra were de- 
tected mostly in hard (probably non-thermal) X-rays in 
the energy range up to 200 keV; the frequencies of this 
QPO range from 18 Hz up to 600 Hz. The QPO emission 
is observed only during certain rotational phases. One of 
several remarkable features of this discovery is the rather 
large amplitude of the luminosity modulation, of the or- 
der of 10% for the high / modes, given that a crust would 
break or melt if the crustal oscillations were too violent. 
The shear strain of crustal oscillations is probably lim- 
ited to about ~ 10%, which, for Z '--^ 10, would result in a 
limit for the oscillation amplitude about one or two per- 
cent of the NS radius Rns, unless vortex pinning and/or 
strong magnetic field make the crust more resilient. 

Even more remarkable is that such large modulation 
appears most prominently even at certain high I har- 
monics. Even if temperature fluctuations of order several 
percent could be somehow achieved at the photosphere, 
it would not lead to such large luminosity modulation if 
many maxima and minima were visible to the observer 
at any given time. Moreover, the high value of I suggests 
that energy in crustal oscillations is triggered from a rela- 
tively small part of the crust's surface and is distributed 
into a large number of modes. It also raises the pos- 
sibility that the deep crustal heating that can account 
for the medium and long term X-ray afterglow follow- 
ing giant flares and other periods of magnetic activity 



(|Lvubarskv et al.ll2002l : iKouveliotou et al.ll2003( ) may be 
powered by crustal oscillations. 

In this paper we propose a model for the amplitude 
variations that appears to support the idea that the oscil- 
lations are indeed torsional oscillations of the crust, and 
may also support the suggestions Thompson ct al. (200^ 
that the magnetospheric currents play a significant role 
in the magnetar emission at high photon energy. 



2. THE MODEL 

We suggest that the modulation of the hard X-ray ra- 
diation visible in the tail of the giant flares is caused by 
modulation of the particle number density in the mag- 
netosphere due to shaking of the magnetic field lines by 
oscillations of the NS crust. The specific candidate emis- 
sion mec hanisms can be cyclotron scattering of thermal 
photons (|Lvutikov k Gavriill [2001 . non-resonant scat- 
tering of thermal photons, or bremsstrahlung off a target 
corona, though a detailed investigation of the latter two 
is not attempted here, as they seem to us less promising. 
The electrical resistance created by photon drag during 
the tail phase of a giant flare is several orders of mag- 
nitude greater - as is the non-QPO photon background 
against which the QPO component must stand out - than 
during the quiet time persistent emission. This consid- 
eration leads us to favor resonant scattering, for which 
the optical depth olfered by the current carriers to the 
thermal photons is the greatest-'^. 

In order to support the twist of the magnetic field, a 
current density j must flow along magnetic field lines. 
The minimum particle density in the magnetosphere of 
the magnetar is therefore j'/(ec). The strength of the 
magnetic field where resonant inverse Compton scatter- 

^ The coronal density may itself be enhanced during tail-phase 
emission by radiative trapping of plasma on closed field lines, and 
this might also be a source of enhanced electrical resistivity that 
directly taps the energy in magnetic field twist. However, here we 
focus only on resonant scattering. 
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ing on electrons occurs is, in the absence of Doppler shift, 



Br, 



0.086£;[t''''lO^'G, 



(1) 



where e'^^^ is the initial energy of the scattered photons 
in units of keV. For magnetars with the surface magnetic 
field 10"'^^ — 10^^ G the resonant surface for hard X-rays 
lies at the distance of the order of ^ 10 stellar radii (i?Ns) 
or less. For the observed frequencies of the QPO the elec- 
tromagnetic wavelength Aqpo = c/vqpo corresponding to 
these oscillations are of the order 50 to 1700 i?NS- Hence 
the resonant surface is well within the wave zone and 
the magnetosphere can be considered as quasistatic. All 
physical quantities changes with time as e~*'^*. 

If the reconfiguration of the m agnetic field induces os- 
cillations of the NS crust (iDunc an 1998), the motion of 
the star surface moves the magnetic field lines frozen into 
the surface. The motion of the field lines induces the elec- 
tric field and, for modes which move the footpoints of a 
given magnetic field line in opposite directions, induces 
a twist in the field line. The twist of the field lines in- 
duces electric current along the lines and hence changes 
the particle number density in the magnetosphere. As 
we show in appendix]^ only modes which move the foot- 
points of a given magnetic field line in opposite directions 
and thus induce the electric current along magnetic field 
lines could modulate hard X-radiation of magnetar. 

The optical depth to the resonant Compton scatter- 
ing on particle with th e charge Ze is (cf. eq. (23) in 
iThompson et all (|2002f )) 



T{e) ^ n'^Zen (1 + cos^ OkB) 



dB 
dr 



(2) 



where all quantities are taken at the corresponding point 
of the resonant surface (rres, ^); nis the resonant particle 
number density, 9kB is the angle between the photon 
wave vector and the magnetic field. For dipolar magnetic 
field \dB/dr\ ~ {3/2)B/r. The current density induced 
by the twist of the magnetic field line is 

c ^ - c SB 1 i B 

47r 47r l\x I l\x irs 

where SB ^ 2{£^/s)B is the variation of the magnetic 
field due to twist of the magnetic field line; Ax is the 
characteristic scale of the electromagnetic field variation; 
s is the length of the corresponding magnetic field line; ^ 
is the oscillation amplitude. Assuming that the plasma 
density is proportional to the current density, see eq. (|Bip 
in appendix [B1 we get for the optical depth 
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where k is multiplicity of pair plasma. The maximum op- 
tical depth is achieved near the equator, where s/r^cs — 
2.6; for oscillation mode with the harmonic number I 
Ax ~ 7ri?Ns/'- Using these r can be estimated as 



kI- 



3i?NS 



K—- 



8 i?NS 



(5) 



Hence the modulation of the optical depth of the order 
of 0.01k would require oscillation amplitude for the har- 
monic number ^ ~ 10 of the order of 1%. Such oscillation 
would produce shear strain of the NS crust of the order 



of 10%, implying that the nearest neighbor distance os- 
cillates by about 1% (In appendix|B]we argue that k may 
be of the order of a few or larger, which would provide 
optical depth modulation of the order of ~ 10%). The 
question is whether such a crustal strain could accommo- 
date the observed QPO in the flux. We now do a more 
quantitative calculation relating optical depth to crustal 
strain and flux modulation. 

3. AXISYMMETRIC FORCE-FREE MAGNETOSPHERE OF 
AN OSCILLATING NEUTRON STAR WITH DIPOLE 
MAGNETIC FIELD 

Consider an idealized case of a non-rotating star with 
dipole magnetic field oscillating in toroidal oscillation 
mode with small amplitude. We consider disturbance 
of the magnetosphere by toroidal modes with even value 
of / and 171 = 0. We use a spherical coordinate system 
centered at the NS center, with z-axis parallel to the 
dipole momentum and the same representation for mag- 
netic field and current s in the force-free magnetosphere 
as in lTimokhi"iil (|2006f ) . For an axisymmetric force- free 
field, the toroidal component induced by the oscillations 
is 

B, = -^,, (6) 
c r sm 

where 27r7 is the total induced poloidal current flowing 
between the pole and the colatitude 9. For small pertur- 
bation of the magnetic field we have 
r sin 9 dcf) B^, 



ds 



B, 
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pol 
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where -Bpoi is the poloidal component of the magnetic 
field and ds is the distance along the field line. The total 
twist of the magnetic field line (the difference — 02 of 
the field line footpoint azimuthal angles) 54>* is 



d9 



d9, 



(8) 



where 9^, is the colatitude of the magnetic field line foot- 
point in the upper hemisphere. The derivative d(j)/d9 can 
be expressed through the derivative d<p/ds as 



d9 



d(j) ds 
'dsM 



d<p sm 9 
ds sm 0, 



VT+3cos2^, (9) 



and d(j)/ds can be obtained from eq.©. Expressing B^ 
through the poloidal current /, eq. and performing 
the integration, we get for the total induced poloidal cur- 
rent 

/=^'50*^Soi?NS^. (10) 

2 Stt cos 9^ 

The total twist of the magnetic field line S(j)^, is 

where ^{6^) is the displacement of the magnetic field line 
footpoint. For the toroidal oscillation modes under con- 
sideration, the displacement of the NS surface 5(6'*) is 

e(^o--eo^^, (12) 

where is the mode amplitude and Y/o is the spheri- 
cal harmonic of order 1,0. The current density can be 
expressed through the magnetic fiux function ijj as 

J = = —.Bpo\ (13) 



r sin0 
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If the resonant surface is located not very far from the 
NS, the optical depth has several maxima at different co- 
latitudes. The position of these maxima as well as their 
amplitude will depend on the shape and position of the 
resonant surface. The latter depends on the magnetic 
field strength as well as on the energy of the upscattered 
photons. In Fig. [T] we plot three resonant surfaces for 
different energies of the upscattered photons for the po- 
lar magnetic field Bq = 10"'^^ G. It is evident that for 
different photon energies the angular distribution of the 
optical depth will be different. In Fig.[2]wc plot the angu- 
lar distribution of the optical depth for different energies 
of the upscattered photons, different oscillation modes 
and different magnetic field strength. In each of these 
plots we assumed the amplitude of the oscillation mode 
^0 = 0.01i?NS- Latitudinal distribution of the oscillation 
amplitude ^ (6', ) as well as the largest coniponen t of the 
strain tensor ee,p (see e.g. iKosevich et"all 119860 for the 
considered modes / = 4; 8; 10 are shown in Figs. [3] and [4] 
correspondingly. 

The optical depth distribution has the following prop- 
erties 



Fig. 1. — Current density distribution in the magnetosphere 
of a NS with dipole magnetic field oscillating in toroidal mode 
(Z = 10, m = 0). The intensity of color in each point is proportional 
to |j(r, 9)1^/^ (this moderates contrast of the plot). Positive values 
of the current density are shown by red colors, negative ones by 
the blue color. By the black lines the resonant surfaces are shown 
for three photon energies Eph = 5; 10; 30 keV - solid, dashed and 
dotted lines correspondingly. The magnetic field at the pole is 
assumed to be Bq = 10^** G. 

The magnetic flux function tp for the dipole field is 



BoR^ 



NS 



dl 



z /• 

where Bq is magnetic field strength at the pole, 
derivative dl/dip is 

dl _ dl d9^ _ 
dip d6^: dip 

Combining equations 

-3 

I f \ 

J = Jo 



dO^ 



BqR^^ sin( 



, cos( 



TSl) we obtain 



where 
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As an example, we plot in Fig. [T] the current density 
distribution for toroidal mode I = 10. The current num- 
ber density is strongly modulated in latitudinal direc- 
tion and therefore the resonant Compton optical depth 
is highly anisotropic. 

Now we can obtain an expression for the optical depth. 
Substituting the expression for the current, eq. 
the formula for the optical depth ([2]) , we get 



mto 



S.0 

'^24i?Ns''''°^l 



1-1-7 cos^ 



3 cos^ 

Here for simplicity, as in [Thompson et al.l (|2002f) . we as- 
sume radially streaming photons. 



• The latitude of peak optical depth is strongly en- 
ergy dependent for the low latitude peaks and less 
so for the peaks near the polar axis (see Fig.[T|). 

• The maximum optical depth does not depend sig- 
nificantly on the initial photon energy (Fig. [21 left 
panels). 

• For moderate 1 (Z < 10) the optical depth produced 
by mode amplitude go varies little with I, although 
the angular dependence is naturally different, see 
Fig. [2] (middle panels). For larger I equation ([5]) 
gives a rather good estimate for the maximum value 
of the optical depth, which, for the same mode am- 
plitude at the magnetic field line footpoint ^(6**), 
grows with increasing I. 

• The maximum value of the optical depth does not 
depend significantly on the strength of the magne- 
tar magnetic field for commonly accepted parame- 
ters of Bq = lOi-* - 10^5 G, see Fig. [21 right panels. 

We consider now three possibilities for modulation of 
hard X-radiation: a) modulation is independent on pho- 
ton energy change and depends only on anisotropic ob- 
scuration of photons generated within the resonant sur- 
face; b) modulation is achieved by upscattering of pho- 
tons in energy by either (i) non-relativistic electrons or 
(ii) at least moderately relativistic electrons. 

Case a) requires an anisotropic optical depth of the or- 
der of 10%. This is easily achieved with the multiplicity 
K of the order of 10, such multiplicity is easily achieved 
because pair-production requires primary electron ener- 
gies of the order of ~ lOmc^, where m is the electron 
mass, (see appendix [B|) . 

For the case b(i) the fluctuation of the photon fiux is 
increased by the Compton-Getting effect as 



U 



d\ogFp^{Ep 
d log -Bp. 



c 



(20) 



where U is the characteristic electron velocity at the scat- 
tering point (the sum of bulk, thermal and drift velocity) 
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1=10; Bo = lxlO'*G Bo = lxlO'*G; Ep|, = 10keV 1=10; Epi, = 10keV 

Ep,,=5keV Eph = 10keV Ep,,=30keV 1=4 1=8 1=10 B|, = lxlo'*G Bo=3xlo'*G Bo=9x10'''g 





Fig. 2. — Angular distribution of the optical depths to resonant Compton scattering for radially moving photons in the magnetosphere of 
an oscillating NS for oscillations with the amplitude go = 0.01-Rns- the upper figures the angular distribution t{9)/k is shown in polar 
coordinates (t{8) / k, 8), in lower figures the same distribution is shown as a plot t{9) for 9 £ [0,90°]. In the left figures the optical depth 
is shown for 3 different initial energies of the scattered photons (-Eph = 5; 10; 30 keV), the magnetic field Bq = 10^* G and the harmonic 
number of the toroidal mode I = 10 are fixed. In the middle figures the optical depth is shown for 3 different modes [l = 4; 6; 10), the 
magnetic field Bq = 10^'^ G and the initial energy of the scattered photons i?ph = 10 keV are fixed. In the right figures the optical depth 
is shown for 3 different values of the NS magnetic field (Bq = 1; 3; 9 X 10^'* G), the harmonic number of the toroidal mode I = 10 and the 
initial energy of the scattered photons i?ph = 10 keV are fixed. 



and Fp\^ is the number density of photons with the en- 
ergy i?pi,. We have used the time integrated photoiisgec- 
trum Fph(-Eph) of the tail phase reported bv iHurlev et all 
([2n05h . This is more than the contribution of the case 
a) if J7/c > 1/4. It gives at least 6% modulation of the 
hard X-ray flux for displacement amplitude of the order 

of 0.01i?NS- 

Case b(ii) minimizes the optical depth implied by the 
twisted field, but, on the other hand, it enables the scat- 
tering process to lift photons directly from the thermal 
peak at E^i, ~ 10 keV to the high energy tail. We es- 
timate that the unmodulated ta il repres ents about 0.1 
of the total photon number fsee iHurlev et al. 2005), so 
a QPO modulation amplitude in the high energy tail of 
10% would require an optical depth to resonant inverse 
Compton scattering of about 1%. 

For oscillation modes with moderate I < 10 the max- 
imum value of t/k reaches ~ 1% (1.5)% when the dis- 
placement of the NS crust is 0.6% (1%), see FigEl The 
strain of the crust caused by a particular mode does not 
exceed 10%, see FiglH Not surprisingly, this is close to 
the maximum beyond which the crust would either break 
or plastically deform, but it seems tolerable, unless there 
are so many modes excited that the root mean square 
strain is much more than per individual mode. The 
strong magnetic field may provide a rubbery consistency 
to the crust material so that it could withstand larger 
reversible strain than brittle crystal, and this would be 
an implication of the likely possibility that many other 
modes were excited in addition to those that produced 



observable QPO behavior. 

Angular modulation of the optical depth to the reso- 
nant Compton scattering may result in temporal modula- 
tion of the hard X-ray photons. Indeed the QPO features 
are ob served only at specific phases of the X-ray l ight 
curve (jStrohmaver fc Watts! [200l [Israel et al.ll2005f) . A 
quantitative fit to the light curves is beyond the scope 
of this article, but we suspect that angular dependence 
of the optical depth is the main reason that the QPO 
features are rotational phase dependent. Altogether, 
the proposed scenario seems to be in good agreement 
with the interpretation of the largest amplitude (~ 10%) 
QPOs as somehow resulting from NS oscillations with 
harmonic number I < 10 (e.g. lStrohmaver fc Wattsir2006l : 
llsrael et al.ll2005| ). 

4. DISCUSSION 

In all of the above, we have assumed for simplicity that 
the torsional oscillations and magnetic field are coaxial 
and axisymmetric. In reality, a more complex topology 
of the magnetic field is likely, hence a more complicated I 
dependence. The general principle nevertheless remains 
that some modes twist field line more than others and 
that some field lines will be twisted more than others for 
any given mode. 

Our proposal that the QPO oscillations result from 
oscillation-induced enhancements of the magnetosphcric 
currents is based on the premise that the emission or up- 
scattering increases with density. This would not work, 
for example, for blackbody emission. We therefore con- 
sidered it significant that these QPO were observed in the 
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Fig. 3. — Displacement amplitude of the NS surface normalized 
to the NS radius as a function of the colatitude for toroidal oscil- 
lation modes Z = 4 — solid line; 1 = 8 — dashed line; / = 10 - dotted 
line. The mode amplitude Jo = 0.01-Rns- 
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Fig. 4. — The largest component of the strain tensor, 2eg^, as 
a function of the colatitude for toroidal oscillation modes I = A - 
solid line; Z = 8 - dashed line; Z = 10 - dotted line. The mode 
amplitude Jo = O.OliJt-js- 

hard X-ray band and not in the thermal peak. Our model 
predicts that QPO behavior should be significantly less 
in any black body thermal component. Similarly the fact 
that QPO are seen at high photon energies but not at 
low photon energies favors upscattering over anisotropic 
obscuration. Our interpretation of the QPO peaks in 
the tail phase emission of the giant flares seems consis- 
tent with the overall picture of magnetar emission that 
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has emerged over the past several years. In particular it 
is consistent with the view that the hard X-radiation is 
due to non-thermal reprocessing by the magnetospheric 
current carriers or direct emission from a hot corona. It 
is also consistent with the interpretation that these os- 
cillations are torsional rather than some other kind. 

One unsurprising implication of the remarkably large 
amplitude of the QPO modulation is that the crust re- 
ceives an enormous amount of mechanical energy from 
the giant flare. We estimate over 10^^ erg per mode if 
the amplitude of oscillation is 10"^ neutron star radii and 
the QPO frequency centered on 150 Hz. Assuming that 
the the modes are indeed torsional modes, it follows that 
their mechanical energy deposition extends throughout 
the entire crust. The fact that the measured Q of these 
oscillations is always much less than the number of pe- 
riods that the QPO behavior lasts is quite curious, and 
the most obvious interpretations are either that a) the 
oscillations damp on a short time compared to the over- 
all duration, requiring continuous excitation, b) many 
modes are excited {'3> 10, as at least 10 separate fre- 
quencies are detected in the individual source) with an 
unresolved spreads in frequencies clustering around the 
observed central ones. Either interpretation implies that 
far more than 10^^ ergs are excited, we may suppose that 
over 10^^ erg is released in the inner crust as mechani- 
cal motion. This is far more than the emission during 
the tail phase and probably goes mostly into heating the 
crust. This is consistent with the assumption that over 
10^^ erg can be d eposited in the deep crust that was made 
(|Kouveliotou e t al. 2003; Eichler et al. 2006) to account 
for observed long term X-ray afterglow from magnetars 
1627-41 and 1900-1-14. 

In each case a decent fit to the observed cooling 
curve requires a large enough mass of the NS that the 
direct Urea process cools the core (at least 1.4Mq). 
Combining this resul t with analysis shown in Fig. 5 in 
IStrohmaver fc Watts! (|2005D would set strong limits on 
the equation of state of the NS. 
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APPENDIX 

A. GOLDREICH-JULIAN CHARGE DENSITY 

The twist of magnetic field lines induces currents flowing along the field lines. Additionally, motion of the field lines 
induces an electric field, and, in order to keep the magnctosphere force-free the charge density in the magnetosphere 
must contain an oscillating component - the " oscillational" Goldreich- Julian (GJ) charge density. Let us compare 
minimal particle number densities required by changes in the GJ charge density and by current supporting twist of 
the magnetic field lines. 

The oscillational GJ charge density could be estimated as 

PG. = ^V-i?^— -— (Al) 

47r 47r c i\x I i\x Aq^o 

Here Ax is the characteristic scale of the electromagnetic field variation, v is the amplitude of the oscillational velocity, 
^ is the oscillation amplitude. The current density induced by the twist of the magnetic field line is given by eq. ©. 
The ratio of the corresponding particle number densities is thus 

^twist Jtwist 1 _ '^QPO ^ -j^ (A2) 



the last inequality follows from estimates of the sizes of the resonant surface and the wave zone given in section [2l We 
note that this estimate is valid for any point in the magnetosphere, also for a point at the resonant surface, s is then 
the length of the magnetic field line from the considered point to the corresponding symmetrical point at the resonant 
surface. Taking into account that the resonant surface is at distances of the order of ~ 4 — f Oi?NS we conclude that the 
particle number density modulation due to changes in the GJ charge density is always less than the modulation caused 
by twist of the magnetic field lines, even for oscillation with frequencies of the order of 1 kHz, when the wavezone size 
is 30i?NS- 

The magnetic field lines get twisted only if their footpoints are moving in opposite directions. Hence, in an axisym- 
metric background magnetic field only half the oscillation modes twist the field lines. The current induced by the twist 
flows along the field line and changes the particle density in proportion to the local value of the magnetic field, so 
the current density with the distance from the NS would decrease as ^ r^^ . These are the modes, we propose, which 
produce observable quasiperiodic modifications in the hard X-radiation in the tail of the giant flares. 

The modes with different symmetry, when the field line footpoints move in the same directions produces only 
modulation of the particle mimber density due to changes in the GJ charge density. Such situation was considered in 
[Timokhin et al. (2000); Timokhin (2007), where analytical 3D model of the force- free magnetosphere of an oscillating 
NS was constructed under assumption that the current along magnetic field lines induced by the stellar oscillations is 
negligibly small. In that case only the charge density near the NS is modulated, which as we shown above introduces 
much smaller modulation to the particle number density comparing to the modes twisting magnetic field lines with 
the same oscillation amplitude. Moreover, far from the star, near the resonant surface, the alteration of the particle 
number density is very weak because pcj falls very rapidly with the distance from the NS, as for a mode with 

harmonic number I. The changes in the particle number density near the resonant surface, in that case, is negligibly 
small and such modes will not produce detectable modulation of the hard X-ray emission. 

B. PLASMA PRODUCTION IN THE MAGNETOSPHERE 

At the tail phase of the giant burst, a magnetically trapped fireball radiates the energy stored during the initial 
outburst. Luminosity at this stage, ~ 10^^ erg/s, is highly super-Eddington. More exactly, the fireball emits radiation 
in the extraordinary mode, for which radiation opacities are st rongly suppressed by magn etic field, and the fireball 
luminosity is in fact the Eddington luminosity for this radiation ([Thompson fc Duncanl[l995| ). However, extraordinary 
photons with the energy E^^ > 40 keV are split in the magnetosphere to ordinary mode photons, for which the 
magnetized opacities are comparable with the non-magnetized ones. Therefore any plasma in the magnetosphere will 
be pushed upwards by the radiation from the fireball. However, the magnetosphere could not become empty because it 
should support electric currents generated during the initial outburst when the the magnetosphere was not force-free. If 
the plasma density falls below the value necessary to maintain the cu rrent, the electric field (displacement current) arises 
and a ccelerates particles until they produce electron-positron pairs (iThompson et al.]|2002l : iBeloborodov &: Thompso"nl 
[2003). 

During the tail phase, the pairs are easily produced via 7 — 7 conversion of photons upscattered off relativistic 
electrons. The photon density in the magnetosphere is n^^ — I0^^L4i/£5o cm~^, where L — 10^^L4i erg/s is the 
luminosity at the photon energy e = 50£5o keV. The scattering cross-section of the ordinary photon is ct = ctt sin^ 6' , 
where 6' is the angle between the photon and the magnetic field in the proper electron frame. So the electron 
with the Lorentz factor 7 = lO/e^o produces a photon with the energy j'^e — lOmc^/e^o after passing a length 
/ = (an)^^ = 20e5o/i4i cm. This photon is easily converted into a pair by collision with a background photon. In 
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principle, electrons with lesser Lorentz factors could also produce pair by upscattering photons with higher energies 
however, the density of photons sharply decreases with the energy therefore one can expect that the electron is 
accelerated to at least 7 > 10 before it produce a pair. After the first pairs are produced, they screen the electric field 
and acceleration ceases. The accelerated electrons rapidly loose their energy upscattering photons. As any ordinary 
photon with the energy > 1 MeV is converted into a pair when it acquires an appropriate angle with the magnetic 
field, one can expect that eventually any primary electron produces at least a few pairs, the number of pairs being 
proportional to the energy the electron acquires within the gap. An important point is that the number density of the 
primary electrons is proportional to the current density, npiimary = j/ec, because there are no other particles in the 
gap. Therefore one can expect that the total number of the pairs is proportional to the current 

n = Kj / ec] (Bl) 

where the multiplicity k is at least of the order of a few but may be larger. In order to find the multiplicity one 
should solve self-co nsistcntly electrodynamics of the gap together with the particle dynamics and pair-production 
(jBeloborodov fc Thompson. 2007,) . which is beyond the scope of this letter 



